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We provide an analysis of basic quantum information processing protocols under the effect of 
intrinsic non-idealities in cluster states. These non-idealities are based on the introduction of ran- 
domness in the entangling steps that create the cluster state and are motivated by the unavoidable 
imperfections faced in creating entanglement using condensed-matter systems. Aided by the use of 
an alternative and very efficient method to construct cluster state configurations, which relies on 
the concatenation of fundamental cluster structures, we address quantum state transfer and various 
fundamental gate simulations through noisy cluster states. We find that a winning strategy to limit 
the effects of noise, is the management of small clusters processed via just a few measurements. 
Our study also reinforces recent ideas related to the optical implementation of a one-way quantum 
computer. 
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I. INTRODUCTION 

The one-way quantum computer (QC C ) [1, 2] is an 
intriguing model for quantum computation (QC) which 
represents an alternative viewpoint from the canonical 
quantum circuit model [3]. The idea behind QC C schemes 
is that single-qubit rotations and two-qubit entangling 
operations can be effectively simulated by the use of an 
entangled resource, the so-called cluster state [4] shared 
by a group of qubits. Quantum Information Process- 
ing (QIP) is carried out by an adaptive set of single- 
qubit measurements, specifically designed for the par- 
ticular simulation desired. Until present, several stud- 
ies have been undertaken in relation to the generation, 
characterization and detection of cluster states, as well 
as more general graph-states [4-7]. In addition, it has 
been recognized by Nielsen [8] that cluster states could be 
used in order to improve the efficiency of the scheme for 
linear optics QC suggested by Knill, Laflamme and Mil- 
burn [9] . Very recently, we have seen the first experimen- 
tal evidence of one-way computation, with the realization 
of a quantum search algorithm achieved through an all- 
optical setup [10]. The class of cluster states and their 
related protocols for QC most certainly deserve a deeper 
analysis and understanding. For example, it is still not 
clear if the one-way model is more robust against the ef- 
fects of external sources of noise and imperfection than 
the quantum circuit model. Several attempts with re- 
gards to this have shown that any Markovian noise in the 
one-way model can be mapped onto non-Markovian im- 
perfections affecting the corresponding quantum circuit 
being simulated [11]. Thus, one can take advantage of 
already existing theorems, valid for the quantum circuit 
model [12], to find proper thresholds for fault-tolerant 
cluster state computation [11]. However we have recently 
introduced an experimentally motivated model for imper- 
fect cluster state generation, which paves the way toward 
the study of the performances of one-way computation 



exposed to realistic sources of non-ideality [13]. Several 
important results were presented, which showed how the 
effectiveness of simulations are affected by randomness 
introduced into the structure of the cluster state being 
used. This led us to search for schemes of cluster state 
QC based on the use of smaller cluster configurations 
than originally devised [1]. In this paper we analyze 
the effects of the model for noise on various key QIP 
protocols, while at the same time using a concatenation 
technique for the construction of compact configurations 
for cluster state QC. Elementary cluster basic building 
blocks (BBLVs) will be introduced, such that they can be 
concatenated to form larger structures for specific circuit 
simulations. With this alternative technique, one may 
bypass the cumbersome stabilizer formalism approach [1] 
and quickly construct more economical circuit simula- 
tions using cluster states. 

The paper is structured as follows: In Section II, we re- 
view the basic properties of cluster states and introduce 
the notation which will be used throughout the rest of 
the paper. The concatenation approach will be discussed 
in Section III and compared to the standard approach 
of constructing cluster state configurations. This alter- 
native viewpoint will then be adopted throughout the 
rest of our analysis. Section IV introduces our model 
for imperfect cluster state generation, where we review 
the fundamental properties of noisy cluster states found 
in [13]. The implications of the profound differences be- 
tween the underlying ideal and non-ideal cluster state 
structures, supplies us with strong motivations to seek 
a deeper analysis of QIP protocols using noisy cluster 
states. Section V is devoted to the transfer of quantum 
information through a linear cluster of arbitrary length 
and we see how the average transfer fidelity depends on 
the dimension of the cluster, with this rapidly becoming 
smaller than the best classical fidelity. Such a result is 
seen as an effect of spreading pre-existing noise by the in- 
creasing number of single-qubit measurements required 
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in order to process quantum information encoded in a 
cluster state. We call this phenomenon inheritance of 
noise by the surviving qubits in a cluster after a given 
set of measurements. This effect is further investigated 
in Section VI, where we address the simulation of singlc- 
qubit and two-qubit entangling quantum gates in noisy 
cluster states. Our analysis includes the basic gates, such 
as the controlled- NOT (CNOT) [3] and general single- 
qubit rotations. The results of our investigation lead to 
the conclusion that QIP in the presence of noise is viable 
only for registers of just a few qubits. This pushes us to- 
ward research of more economical strategies for gate sim- 
ulation [13] than those suggested in the original one-way 
model [1]. Finally, Section VII summarizes our results. 



II. BACKGROUND FOR CLUSTER STATES 

We begin our investigation by reviewing some of the 
elements of the QC C model. The stabilizer formalism will 
be used extensively in this Section. This is a necessity, as 
we will want to compare it with the more advantageous 
cluster state concatenation approach. 

A cluster state \4>{k,}) c is a pure multipartite entan- 
gled state of qubits positioned at specific sites of a lattice 
structure known as the cluster C. The complete charac- 
terization of cluster states is the starting point in under- 
standing of how the QC C model works. A cluster state 
\^{k}) c is formally defined as the eigenstate of the set of 
Hermitian operators 

#<«*>= <#> (g) o-W (l) 

6Gnghb(o)nC 

such that K^\4> {K} ) C = (-1) K '\4> {K} ) C [1, 2]. Here, a x 
(a z ) is the x (z) Pauli matrix, while each acts on 
the qubit occupying site a G C and on any other qubit oc- 
cupying a neighboring lattice site b G nghb(a) fl C. They 
are correlation operators which form a complete set of \C\ 
independent and commuting observablcs for the qubits 
within C. Thus, a cluster state |0{ K }) C is completely spec- 
ified by the set of numbers {k} := {n a G {0, 1}| a G C}. 
For example, a different cluster state \4>{k'}) c belonging 
to the same physical cluster C (i.e. the same lattice con- 
figuration), is completely described by the elements in 
the set {«'} := {n' a G {0, 1}| a G C}. All the clus- 
ter states of an N qubit cluster Cm correspond to the 
2' c l possible combinations of the elements in {k,}. They 
are mutually orthogonal and form a basis in the 2 N - 
dimensional Hilbert space of the cluster. These states 
are also all equivalent under the application of af^ on 
individual qubits. In this paper, we will often use the 
simplified notation \4>)c '■= |0{k})c; with {k} = {0, Va G 
C}. This particular cluster state is generated by first 
preparing the product state \+)c = ® a ec °^ * ne 
qubits at all the sites a of C. In our notation, |±) a = 

(1/V2)(|0) ± |l)) a are the eigenstates of a ( x a) and |0) o 



(|l) a ) is the eigenstate of a z a ' corresponding to the +1 
(— 1) eigenvalue. We then apply the unitary transforma- 
tion — ]\( a b ^ S ab to this initial state, where (a, b) := 

{a,beC\b-ae ld } and 7l = {1}, 7 2 - {(1, 0) T , (0, 1) T }, 
73 = {(1, 0, 0) T , (0, 1, 0) T , (0, 0, 1) T } for the respective di- 
mension d of the cluster being used. Each operator S ab 
can be described by the controlled phase gate 

S ab = |0) a (0|®l (b) + |l) <l|®trW (2) 

and is an entangling operation between the qubits at sites 
a and b of the cluster C. It is important to note that all 
S ab, s mutually commute and therefore the time required 
for the generation of a ./V-qubit cluster state is indepen- 
dent from N [1]. The state generated by the action of 
S< c ) on \+) c is found to be [1] 

S^\+) C ^ 1] ^ h (g)|+)aH</>}c. (3) 

<a,b> aeC 

An Ising type interaction between the qubits in a square 
lattice produces a time-evolution operator having the 
form of S( c \ 

A. QC'c Operation 

Let us assume we want to carry out a unitary operation 
U g associated with a particular quantum gate g which 
acts on an unknown input state \ipi n ) of n logical qubits. 
The idea behind the QC C model is to use a cluster of 
N qubits in a particular physical configuration, C(g). In 
order to better understand how the U g operation on C(g) 
takes place, it is convenient to consider the cluster as 
partitioned into three sections. First, we consider an n- 
qubit input section Ci(g) which encodes the input state 
|V>in) of the register. We then have a body section Cm (<?) 
and an n-qubit output section Co(g) for the read-out 
of Ug. The three sections have no mutual overlap and 
taken altogether, they reconstruct the structure of the 
cluster (i.e. C a (g) nCp(g) = 0, (J a C a (g) = C(g), where 
a,p € {I,M,0} and a ^ (3). We begin by preparing 
iV'in) on the cluster qubits in Ci(g) and denote this state 
as \ipin)ci(g)- ^ satisfies 

n 

IV>i„)eH 9 )=^>*](g)l+}, (4) 

with {Z} := {zi\i G 2™} and the multi-qubit projector 
Pz'lzy\ a i] = 2i=i a iPz'^ ■ Here, we have introduced 
the single z, state projector 

which projects each qubit in Ci(g) into the state \s^)j, 
with as the value of the j-th binary digit of the inte- 
ger Zi and [j] corresponding to an operation on the log- 
ical qubits of the cluster. It should be made clear that 
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the multi-qubit projector P% does not correspond 

to measurements which are performed in reality, but is 
fictitiously introduced in order to relate the logical in- 
put state IV'in) to any possible state present on the input 
qubits of the cluster [1]. In other words, Ci{g) is anal- 
ogous to the register in the quantum circuit model and 
\ipin}ci(g) i s an arbitrary superposition of computational 

states IV'in)cr(s) = Eili"il z i). (EiliW 2 = !)■ Af- 
ter this preparation of \i>in) Ci(g): tne cluster is entangled 
by the operation S c ^ given in cq. (3). A measurement 
pattern M Cm ^ is then applied to the body section of the 
cluster Cm(sO- The pattern is specified by a set of vec- 
tors f a (a € Cm (g)) defining the bases in the Bloch sphere 
of a set of single-qubit measurements to be performed in 
Cm{9)- We denote {s} — {s a e {0, 1} | a e Cm(s)} as the 
set of measurement outcomes obtained after M Cm ^ is 
applied. Performing a measurement pattern on the body 
section of a cluster is formally equivalent to applying the 
projector 

pC^ {M)= l + (-l)^.**> (6) 
keC M (g) 

to the initial state. It is easy to see that S c ^ and 
Pz { ^\cii\ commute. With P^ b) (M) acting only on 
the body section Cm (<?) of the cluster, the state of the 
entire cluster C(g) after this procedure can be written as 

Pz'^Wi) P C { :f 9 \M) S (c ^\+) c(g) . (7) 



While the body qubits have been removed by the set of 
measurements in J\4 Cm ^ 9 \ Eq. (7) still contains degrees 
of freedom in the input section. To eliminate them and 
obtain just the cluster qubits of Co(g), we perform mea- 
surements in the oveigenbasis on each qubit in Ci(g). 
That is 

IVw) Co(9) = Pl£\x)P^ } [amc {g y (8) 

A theorem central to the cluster state model [1] states 
that if \ip)c(g) satisfies the set of eigenvalues equations 

oZ'WJ {U g a^uff° {9) M c(fl) = (-l) A -'|V>c (fl ), 
a c z '^ (y g a^uf) C ° {9] V> C(fl) = (-l) X ^mc(g), 

(9) 

the input state IV'in) and the output |VWt) are related 
via 

IVout) = UgUv IVin), (10) 



where is a byproduct operator, acting locally on the 
logical qubits, given by 

U* = ®(aWy*+^(o-M)^. (11) 

Thus, the realization of the appropriate measurement 
pattern effectively simulates the action of the gate g. 
However, we must counteract the unwanted effect of the 
byproduct operator which appears in the simulation. If 
g belongs to the Clifford group [14], the propagation of 
Us through U g in Eq. (10) leaves U g unaffected [1]. That 
is: U g Uj: = Uj;U g with Uy, = UgU-^Ug 1 which remains 
a product of local operations on the qubits in Co{g)- 
We then just apply the decoding operator to the out- 
put cluster qubits in order to wash out the effect of the 
byproduct operator. On the other hand, if g is not in 
the Clifford group, we obtain U g U^ = UsU g , while still 
keeping the local nature of the byproduct operator. In 
this case, the measurement pattern must be changed in 
order to account for the propagation effect [1]. 

III. CONCATENATION OF BASIC BUILDING 
BLOCKS 

It is clear that the important elements in the 
QC C model are the design of the appropriate measure- 
ment pattern, along with a complete knowledge of the 
corresponding decoding operator. Within the stabilizer 
formalism, the pattern JvI Cm ^ suitable for a gate g is en- 
gineered by manipulating the eigenvalue equations which 
characterize the cluster state [1]. For large N and n, it is 
not hard to believe that this procedure becomes rapidly 
unmanageable and an alternative approach to the con- 
struction of exploitable cluster configurations is required. 

This necessity is justified under practical perspectives 
as well. Indeed, as we discuss in the next sections, the 
realization of a given cluster state is a system-dependent 
issue. A method particularly suited for condensed-matter 
systems (as well as optical lattices, discussed in Sec- 
tion IV), is the creation of regular cluster configurations 
from which redundant qubits, not necessary for the com- 
putational steps performed and thus superfluous in the 
corresponding cluster's layout, are removed through ap- 
propriate single-qubit measurements. This issue will be 
reprised in more detail in Section VI. In ideal conditions, 
the removals are harmless for the performance of a par- 
ticular gate simulation. But when (intrinsic or external) 
decoherence is introduced in the cluster, this is no longer 
true. One of the key results in [13], is that the number 
of measurements performed on the cluster qubits (either 
within M Cm ^ or instrumental to the cluster state gen- 
eration) should be as small as possible. Thus, tailoring 
a cluster state by cleansing it from the redundant qubits 
may not be an efficient strategy. 

On the other hand, very recent experimental efforts 
in the all-optical scenario have demonstrated that small 
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clusters of just a few qubits can be efficiently produced 
within state of the art technology [10, 15]. In these cases, 
the setup is built so that exactly the layout needed is pro- 
duced and no qubit has to be thrown away. Small op- 
tical clusters can then be mutually connected by gluing 
them as suggested in [16]. Here, we describe an approach 
to cluster state information processing which uses fun- 
damental configurations or BBB's. By concatenating 
them, i.e. placing the blocks next to each other such 
that they overlap, any desired layout corresponding to a 
precise computational task can be constructed and the 
appropriate decoding operator is naturally retrieved. 
This concatenation technique gives us the possibility of 
designing optimal elementary steps, each corresponding 
to one stage of an experiment, which limit the effects 
of the intrinsic model for imperfection introduced in [13] 
and discussed in more detail here. 



A. Basic building blocks and concatenation 

We start by introducing the BBB's and give their 
equivalent network circuits. Our diagrammatic notation 
is such that each qubit will be represented by the symbol 
of a square. The angle a inside the i-th qubit symbol 
identifies the basis Bi(a) := {\+) a , |— }"} in which that 
qubit is measured. Here, |±) a = (l/y/2)(\0) ±e M |l)) 
with a the angle from the +ve x axis in the xy plane of 
the Bloch sphere and s^°^ g {0, 1} is the corresponding 
measurement outcome. Empty squares stand for unmea- 
sured (output) qubits. In the configurations we show 
here, the orientation is irrelevant so long as the topology 
does not change. 

The smallest cluster state we can conceive consists of 
two qubits and can be used to simulate a unitary op- 
eration on one logical qubit \Qi) encoded on a physi- 
cal cluster qubit (thus |Cf(.g)| = 1), as shown in Fig. 1 
(a). This configuration will be denoted from now on 
as BBB\. The overall operation simulated by BBB\, 



(a) 

MEEK) 



(b) 



\Qi) - R~ a - H - \Q[) 



FIG. 1: (a): The BBB 1 layout, (b): The operation sim- 
ulated when qubit 1 is measured in the B\(a) basis and 



xy(a) 



= is obtained. 



when the measurement on qubit 1 gives s\ yi *°^ = as 
the outcome, is shown in Fig. 1 (b), where R~ a repre- 
sents a rotation around the z axis by an angle —a and 
H denotes a Hadamard gate. Due to the probabilistic 
nature of the simulation, it is necessary to apply a de- 
coding operator U^{s^^) — Ux 1 to qubit 2. Using 
the same two-qubit cluster layout, with two encoded log- 
ical qubits \Q\) and \Q^) («-e. here Ci(g) = Co(g), see 
Fig. 2 (a)), we simply simulate a controlled- 7r phase gate 



(CrPHASE) [3], as indicated in Fig. 2 (b). We denote 
this configuration as BBB2. 



(a) 
1 



(b) 



IOi> 



IQ2 



\Qi) 



FIG. 2: The CRPHASE gate simulated by BBB 2 on two log- 
ical qubits \Qi) and \Qi) . 

Our set of BBB's is completed by the introduction of 
BBB3 shown in Fig. 3 (a), where qubits 1 and 3 en- 
code the input states and a measurement is performed 
on qubit 2 in the B2(a) basis. This pattern simulates 
the transformation 



Tbbb 3 — 
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r a (12) 

where the top (bottom) sign corresponds to s% — 

(sg — !)■ P° r a = {Oi 71 "} * e - a measurement 
of the bridging qubit in the ±a x eigenbasis, one ob- 
tains a non-unitary gate. We are particularly inter- 
ested in what is found for a = \, which corresponds 
to an operation decomposed as shown in Fig. 3 (b). 
The decoding operator for this specific gate simulation 
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is Uj, — —i(Jz (8 <j\ 



\Qi) 
2 



(a) 
1 

a 



(b) 
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IQi) 

|Q2> 



FIG. 3: (a): The BBB 3 configuration, (b): The equivalent 
quantum circuit corresponding to the operation simulated by 

BBBj, with a = ty/2 and ] = 0. 

Using this set of BBB's, we can construct more com- 
plicated configurations and in order to see this more 
clearly, we use the following observation from [1]: 

Observation 1. Consider a quantum circuit g associ- 
ated with the unitary operator U g simulated on a cluster 
C(g). Let g be comprised of two consecutive circuits g\ 
and gi on subclusters C(gi) and C(gi) respectively, i.e. 
9 = gig\ and C(g) = C(gi) U C(g 2 ), with C(gi) H C(g 2 ) 
containing one cluster qubit for each logical qubit. These 
subgates have associated unitary operators U gi and U 92 . 
The method of entangling the whole cluster C(g) and per- 
forming the required measurements for simulating g, is 
equivalent to entangling the qubits of C{g\), performing 
the required measurements for g\, then entangling the 
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qubits ofC(g 2 ) and performing the required measurements 
for g 2 . 

Proof. In the first method, the logical input state \ipi n ) 
is encoded onto the cluster qubits in the input section 
Ci(gi) of the first subcluster. The entire cluster C(g) is 
then entangled and qubits in C(g) \ Co (32) are measured 
in the required bases, with the logical output state IV'out) 
present on the cluster qubits of Co (32 )■ 

In the second method, the logical input state (Vin) is 
again encoded onto the cluster qubits in the input section 
Ci(gi) of the first subcluster. The subcluster C(gi) is then 
entangled and qubits in C(g\)\Co{gi) are measured in the 
required bases, with the intermediate logical output state 
IV'out) present on the cluster qubits of Co(gi) = Ci{g 2 ). A 
similar procedure for the second subcircuit g 2 is carried 
out and the logical output state IV'out) ^ s present on the 
cluster qubits of Co (32)- 

To show the two methods are equivalent, is to show 
that l^out) = IV'out)- Let Pi, P 2 be projectors repre- 
senting the measurements on qubits in C(gi) \ Co(gi) 
and C(g 2 ) \ C {g2) respectively with S^ c( - 91 ^ =: Si and 
S^-)) = : S 2 representing the entanglement operations 
on subclusters C(gi) and C(g 2 ) respectively. As Pi com- 
mutes with S 2 , we find that 



with R z a a x a and R z a o-y b ■ We may write 



P 2 S 2 PiSi — P 2 PiS 2 Si, 



(13) 



This mixture of entangling operators and projectors acts 
on the state \i>m) Cl ( g ) ®a£C(g)\c,(g) l+) a - As the y act 
equivalently, we have that IV'out) = IV'out) = U g Us iVin) 
and the two methods are equivalent. □ 



Observation 1 can be extended to an arbitrary number 
of sub-gates gi associated with unitary operators U gi on 
sub-clusters C(gi). These then make up the gate g asso- 
ciated with the operator U g on the entire cluster C(g). 
For each projector Pj there is an associated unitary op- 
erator U 9i applied to the logical qubits. However U 3i 
is not the only operation applied. Due to the fact that 
there are two possible outcomes for each qubit measured, 
the gate simulation is probabilistic. It turns out that 
for all the BBB's in the last section, the actual uni- 
tary operation implemented can be written as U gi U-s(si) 
(gi = BBBi, i = 1,2,3), where Ux(si) is a local byprod- 
uct operator dependent on the outcome Si :— §J l,( - Q ^ 
of the measurement performed on the j-th qubit in the 
C(gi)\Co(gi) part of the cluster. In this way the unitary 
operation U 9i becomes deterministic, as we can always 
propagate Us(si) through it, keeping the byproduct op- 
erator local. For a large concatenation of subclusters, 
this propagation procedure can be carried out for U 9i 
and U^(s a ,s b ,s c ) = a s x a a y b a s z c , where s a ,s b ,s c € {0,1} 
are a result of the byproduct operators corresponding to 
the gates acting before U 9i . However, when U 3i is not 
in the Clifford group, it will change upon propagation 
of the byproduct operators through it. For all of the 
BBB's mentioned before, we find that this only occurs 



L jD — a 



~o _ a s bJl-a 



(14) 



where a 1 — (— l) Sa a and a" = (— l) Sb a. In other words 
the measurements become adaptive, but the resources 
(in terms of the number of cluster qubits) never change, 
the byproduct operators always remain local and we 
can counteract the effect of the byproduct operators by 
changing the measurement angle from a to —a. Overall 
in a large concatenated circuit of BBB's, we end up with 
the logical input \ipin) and output IV'out) of the unitary 
simulation U g — Y\a=\n\ U gi , related via 



IV'o 



n 

i=W\ 



n 



Ug, IVin). (15) 



Here, Af is the number of subgates, and the ~ repre- 
sents the fact that the byproduct operators change on 
propagation, or the unitary operations of the subclus- 
ters become adaptive. A final note to add is that some 
byproduct operators will have to propagate further than 
others in this scheme. By applying the hermitian conju- 
gate Ul(s a7 s b ,s c ) = (Ili=|Af| Uz(si)Y of the propagated 
byproduct operators we will recover the unitary opera- 
tion desired. The U^(s a , s b , s c ) are decoding operators 
and for each logical qubit, they will always be a mul- 
tiple from the set {1, a x , a z , o~ y } up to a global phase 
factor € {1,-1}. We therefore recover the cluster state 
model for QC, without the need for the stabilizer for- 
malism and eigenvalue equations of Eq. (9). It is now 
easier to design the cluster configuration which simulates 
a desired quantum circuit and we are on our way to- 
ward the construction of computationally useful extended 
building blocks (EBB's). By concatenating two BBBi& 
and two BBB 2 s, we produce a bidimensional square 
cluster state, the box cluster, which has been experimen- 
tally realized and used in order to perform a two-qubit 
quantum searching algorithm [10]. The physical layout, 
the measurement pattern and the corresponding equiv- 
alent circuit are shown in Fig. 4. On the other hand, 



a 









Qi) 



1) 



IQa) 



RZ 



H 



\Q'x) 



FIG. 4: The concatenation of two BBBis (Fig. 1) with two 
BBB 2 's (Fig. 2). 

by concatenating two BBBiS and a BBB 2 , we obtain a 
simple four-qubit linear cluster which is particularly in- 
teresting. Indeed, if the single-qubit measurements are 
performed in the cr^-eigenbasis, the corresponding equiv- 
alent quantum circuit is locally equivalent to a CNOT 
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gate i.e. the read-out bases of qubits 2 and 4 are the 
a x eigenbases [13]. The layout and the measurement 
pattern are shown in Fig. 5. We will comment on the 



iq0 3 m 



X 


X 


Qi) 2 





\Qi 



\Q» 



h -r 



\Qi, 



H - \Q' 2 ) 



FIG. 5: The concatenation of two BBBi's and one BBB2 
with measurements in the -Bi(O) and -83(0) basis. 

importance of this cluster configuration in Section VI, 
where its role in an economical scheme for cluster state 
based QC is highlighted [13]. For our purposes here, it is 
instructive to explicitly compute the byproduct operator 
corresponding to this gate simulation. As we have said, 
the decoding operator needed after the measurement pat- 
tern in BBB\ is given by either 1 or a x , depending on 
the outcome of the measurement performed. Thus, in- 
cluding the random action of the byproduct operators, 
the transformation in Fig. 5 can be represented overall 

by 



(2) 



hhCrPHASE a 



r (2K 



w(0) 



1< 4 ) H 



(16) 

Note that no rotation operator appears in the above ex- 
pression because of the choice of -Bi(O) and i?3(0) for 
the measurement pattern. We now use the fact that 
C 7r PHASE(o-i i) <g> = (a { x l) cr^QPHASE and that 



H, 



Hi, so that the {a 



y(o) 



1^) part 

in Eq. (16) can be propagated until we end up with 
f7 E H 4 S' 24 H 2 , where 



fr ( 2 ) s r 
Us = a x 



r (4)[»;» 



(17) 



with © as the logical XOR operation in this case only. 
This result, originally presented in [13], is now rigor- 
ously demonstrated by a simple argument based on the 
concatenation technique. Our EBB set is completed 
by the configuration in Fig. 6 (a), i.e. the concatena- 
tion of one BBB3 and two BBB\S. The bases for the 
measurement pattern are B\^(a — 7r/2) and qubits 1 
and 3 encode the input states. The overall transforma- 
tion is decomposed as shown in Fig. 6 (b). Following 
the lines depicted above, the propagated byproduct op- 
erator corresponding to this gate simulation is Us = 

a x 2 <g> g x 2 . We will see in Sec- 

tion VI that this EBB plays a central role in the original 
CNOT simulation discussed in [1]. 

It is remarkable that the explicit form of the byprod- 
uct operators has been found by using the expressions 
valid for the BBB's, the concatenation rule and basic 
commutation properties. The stabilizer formalism is not 
required anymore and the advantages of this alternative 
technique become clear once complicated cluster config- 
urations, involving many qubits, need to be considered. 
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FIG. 6: (a): The configuration obtained by concatenating two 
BBB\s with a BBB3. The single-qubit measurements are 
in the Bi_,3(7r/2) bases, (b): The corresponding equivalent 
quantum circuit. The boxed part is equivalent to a PHASE. 



We thus drop the stabilizer formalism and adopt this 
technique for the remainder of the paper. 



IV. IMPERFECT GENERATION OF A 
CLUSTER STATE 

The key element in generating a cluster state is the 
ability to perform the controlled gates S ab on the qubits 
occupying the sites of a cluster C. However, these oper- 
ations can be inherently imperfect or not precisely con- 
trollable. For example, an optical lattice loaded by a 
bunch of neutral atoms is a candidate for the embodi- 
ment of a QC C [5, 6, 17]. In this particular setup, any 
two-qubit interaction is realized through controlled col- 
lisions of the atoms loading the optical lattice. Several 
experiments [5] consider a unidimcnsional lattice loaded 
by bosonic atoms with two hyperfine levels \c\) and \c2). 
The state of an atom occupying the j'-th site in the lattice, 
prepared in its internal state \c\) (IC2)), can be described 
by the application of the bosonic creation operator cf 

(£2^) to some fiducial initial atomic state. The Hamil- 
tonian describing the atomic interaction takes the form 
of a two-species Bose-Hubbard model [17] 



H i — 



ID* 

3=1 



L 2 /' 



+2c[ 3 "c 



(18) 



where U is the coupling strength [18]. The first two terms 
in the Hamiltonian describe in-site free dynamics of the 
two different species, while the last one takes into account 
the interactions responsible for the controlled phases in- 
troduced in the joint state of two atoms which occupy 
the same site. The coupling strength is in general a func- 
tion of the position in the lattice (i.e. U has a spatial 
profile along the linear lattice) and its value can fluctu- 
ate due to instabilities in the intensities of the lasers used 
in order to create the lattice. A fluctuating U gives rise 
to a random phase shift imposed after the qubit-qubit 
interaction, which may be different from the it value re- 
quired to attain a perfect cluster state. Thus, the effi- 
ciency of the S ab 's critically depends on the control we 
have over the strength of the interactions. The initial 
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filling-fraction of the lattice may also influence the per- 
formances of the gates and the extent to which entan- 
glement can be spread across the cluster [5, 17]. If the 
degree of control is not optimal, we face the problem of 
imperfect (inhomogeneous) interactions throughout the 
physical lattice. Here, this issue is addressed by formally 
considering imperfect entangling operations Sf^'s on the 
qubits of the cluster. Our model consists of controlled 
gates having the form 

S a D b = \0) a (0\ ® l (b) + |1>«<1| ® (|0>6<0| - e^\l) b (l\) , 

(19) 

which add the phase 9 a to the desired and optimal it. 
As in the ideal case, all S£, s mutually commute and 



Il<a,h> 3d 1S 



the imperfect entangling operation Sj-, 
unitary. 

Let us consider the use of Eq. (19) in order to generate 
a unidimensional noisy cluster state (a noisy linear clus- 
ter state). By applying Sp to the product state \+)c of 
N qubits in a linear cluster, we find [13] 



N-l 

\4>)E = S C D \+) C = II (-e ie ¥ z Hzi), (20) 

Zj j = l 

where Zj is the value of the j'-th binary digit of the 
integer z, and the summation runs over all the Zj for 
i € [0,2^ — 1]. The properties of the state \4>)c have 
been discussed in detail in [13], where the fidelity and en- 
tanglement structure of noisy cluster states are analyzed. 
Various probability distributions were placed on the un- 
wanted phases 9j in order to give a better idea of how 
the fidelity might be affected within realistic situations. 
However it was noted that no universal model for the 
probability distribution existed. As a result, we provided 
a more general picture to which specific distributions 
could then be investigated, such as the Gaussian distribu- 
tion. It was found that the fidelity of unidimensional and 
bidimensional clusters rapidly decreases with a rise in the 
number of qubits, even when small amounts of noise are 
considered. The entanglement structure was also found 
to be profoundly affected by the randomness introduced 
to the unwanted phases. It was shown that the amount 
of genuine multipartite entanglement becomes reduced 
to the benefit of bipartite quantum correlations, fragile 
against fluctuations in the unwanted phases 9j and even- 
tually leads to an entanglement breaking effect. From 
this previous investigation, one can clearly see that QIP 
protocols should be strongly affected by imperfections in 
the initial entanglement created between the qubits in a 
cluster [13]. Here we intend to give a broader picture of 
the effect of this noise model on various other key QIP 
protocols. 



Section VI will be dedicated to the basic computational 
steps, here we deal with a communication issue repre- 
sented by the transfer of quantum information through a 
cluster state. We approach the reliability of state trans- 
mission along a quantum channel represented by a linear 
cluster. Quantum state transfer has recently received 
considerable attention in the context of limited-resource 
QIP [21]. Here, our approach will be using a quantum 
resource provided by the multipartite entanglement in a 
cluster state and conditional dynamics given by measure- 
ments. 

The idea behind the realization of information flow 
through a cluster state is simple. In the one-way model, 
seen from the original stabilizer viewpoint, any set of 
single-qubit measurements performed onto the qubits be- 
longing to Cm(<?) ; together with the measurements in the 
cr^-eigenbasis of the qubits in Ci{g) (see Section II A), 
processes the encoded input state and at the same time, 
transfers it to the Co{g) section. Thus, one can naively 
think about simulating the identity gate g = 1 using an 
appropriate cluster state, where the transfer of the input 
state to the output section of the cluster occurs naturally. 
However it is much simpler to think of this information 
flow in a linear cluster in terms of the concatenation of 
many BBBiS. In this case, measuring qubit 1 in the 
f?i(0) basis, means the encoded state is transferred to 2, 
while at the same time rotated by H. Thus, forgetting 
for the moment the effect of the byproduct operator, by 
arranging the concatenation of many BBB^s, i.e. a pat- 
tern of measurements in the single-qubit cr^-cigenbasis 
for all the qubits in Ci{g) U Cm (<?) of a sufficiently long 
linear cluster, either we realize the information flow we 
are looking for (when the number of effective Hadamard 
gates is even) or we simply get |Vw) = H \ip- m ) which 
can be corrected once we know the parity of N . 

The question we raise here is about the efficiency of the 
transfer process when noisy cluster states are used. Ide- 
ally the fidelity for information transfer would be equal 
to unity. However, this is not the case if noisy linear 
clusters are used. In order to obtain a full picture of the 
protocol for information flow, we need to calculate the 
form of the proper decoding operator arising after the 
concatenation of many BBBiS. 



A. Identity and Hadamard gate 



We already know that BBB\ with i?i(0) effectively 

simulates H on a logical qubit |-0i n ), with Us = (Jx 1 ■ If 
we now concatenate two BBB\ 's to obtain a three-qubit 
linear cluster, by measuring qubits 1 and 2 we simulate 



INFORMATION FLOW ACROSS A NOISY 
LINEAR CLUSTER 



ty(0) 



H = <J X 



xy(0) 



11/(0) 



1 



(21) 



(7s 



We now start our analysis of the performances of one- 
way protocols for QIP with noisy cluster states. While 



on the logical qubit \ip{ n ). The decoding operator to de- 
code the information transferred to and stored on qubit 
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FIG. 7: Information flow through an N-qubit linear cluster. 

_ i s *y(0) s xy(0) 

3 is then given by Uj, = a z x a x 2 . The advantage 
of using this concatenation technique with respect to the 
stabilizer formalism is evident and even more striking 
when larger linear clusters are considered. 



B. Flowing information 

In the concatenation model, for an arbitrarily long 
odd-iV qubit quantum channel, the concatenation of 

N — 1 BBBx's results in rii=(iv-i)( cr » < H) applied to the 

input state |'0in), where sf = s* \ Each time the x- 
Pauli matrix corresponding to an even-labelled qubit is 
propagated through an H, thus being transformed into 
a z-Pauli matrix, the product HH = 1 is obtained and 
the x-Pauli matrices of odd-labelled qubits remain un- 
changed. This results in all the even labels (odd labels) 
becoming associated with er z (cr x ) and we obtain the fol- 
lowing structure for the propagated byproduct operator 

fj s f®(Ei=a" 1)/2 4i-i)mod2 sJe(£^- 1)/2 4*)mod2 

(22) 

So far however, we have just dealt with ideal cluster 
states. We now apply the above protocol for t he tran s- 
fer of information encoded in |?/>in) = o|0) + y/l — a 2 \l) 
across an iV-qubit noisy linear cluster. The theoretical 
analysis we have developed so far can still be applied to 
the case of S ab — > i.e. the form of the byproduct 
operator is the same as Eq. (22). We have simulated 
the measurement pattern for information flow on an en- 
coded noisy linear cluster state. Despite the computa- 
tional challenge represented by the exponentially growing 
number of different outcome sets (for an iV-qubit linear 
cluster, there are 2 N different sets {sj}, j € [1, N]), it has 
been possible to explicitly calculate the state transfer fi- 
delity TN(a,9j) = KV'inlVwt)! 2 for each outcome set, up 
to N = 9 qubits. This enables us to evaluate .Fat (a, Oj), 
i.e. the state fidelity averaged over all the different sets 
of measurement outcomes, after the application of the 
relevant decoding operator given by the hermitian con- 
jugate of Eq. (22). This quantity gives us an estimate 
of the average performance of the transfer process. It is 
important to note that this approach is only one of the 
possibilities available, the other being the postselection 
of the event corresponding to a favorable outcome con- 
figuration. For example, one could choose to discard all 
the events but the one where = OVj € [1,N] as this 

case corresponds to = 1^ and therefore no local ad- 



justment is required after the simulation of the identity 
gate. 

However, in order to perform a more quantitative in- 
vestigation and wash out any initial state-dependence, 
the average state transfer fidelity for any input state \ipxa) 
must be considered [21]. This can be done by assum- 
ing a uniform distribution for a and integrating over the 
Bloch sphere surface as Tn{9) = (l/47r) / JF w (a, 0)dS, 
where dS is the surface element. The comparison be- 
tween JFjy(a, 0j) and Fn{9) reveals an almost uniform 
behavior of Fn(o>, 9) with a. On the other hand, the 
integration over the Bloch sphere's surface allows us to 
compare Tn{8) with 2/3, the best fidelity achievable by 
measuring an unknown qubit state along a random di- 
rection and sending the result through a classical chan- 
nel [22] (see also Bose and Paternostro et al. [21]). The 
analysis of TV = 3, 5, 7 and 9-qubit cases in Fig. 8 (a) 
reveals that as soon as 9 ~ 0.65, ,Fg(0.65) < 2/3. By 
increasing 9, all the other transfer fidelities (except the 
case of N — 3) become worse than the classical thresh- 
old value and are thus useless for quantum state transfer. 
The range of exploitability of the cluster channel shown 
in Fig. 8 (a) is considerably small if probability distribu- 
tions attached to the Oj's are considered. We have quan- 
titatively addressed the case of individual Gaussians, cen- 
tered on 9j = and with increasing standard deviation 
a. In Fig. 8 (b), we show the totally averaged fidelity 
Fn,<t « J p(9,cr)TN(9)d9, for a Gaussian p(9,a) which 
retains a parameterisation in N and the standard devia- 
tion a. For a = 1, it can be seen that already for N > 5 
the cluster channel becomes less efficient than the best 
classical strategy for transfer. On the other hand, by 
reducing the amount of randomness in the noisy cluster 
state i.e. by reducing the spread of the distribution down 
to a = 0.5, the usefulness of the transfer protocol is re- 
stored but only for cluster channels of just a few qubits 
(N < 9). Thus, the study of state transfer through lin- 
ear cluster states coincides with the conclusions we have 
drawn about the state fidelity in [13]: whenever intrinsic 
random imperfections are considered within the cluster 
state model, the dimension of the register plays a critical 
role in the efficiency of the QIP protocols. 



VI. COMPUTATION WITH NOISY CLUSTER 
STATES 

We are now in a position to investigate the perfor- 
mance of QC C within the framework of noisy cluster 
states. In this Section, we address both singlc-qubit 
rotations and two-qubit entangling gates, paying par- 
ticular attention to the paradigm of this class of gates, 
the controlled- NOT (CNOT). The gate fidelity relative to 
these operations will be thoroughly studied against the 
effects of the unwanted phases introduced by S^f's. 
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FIG. 8: (a): Fidelity of information transfer with a averaged 
over the single-qubit Bloch sphere and all phases 6j — 9 (Vj). 
From top to bottom curve, we show N = 3 (♦), N = 5 (~k) 
N = 7 (■) and N = 9 (A). The horizontal line (•) repre- 
sents the classical threshold 2/3. (b): Fidelity of information 
transfer with a averaged over the Bloch sphere and 9j aver- 
aged over Gaussian distributions, centered on 6j = 0, with 
standard deviation a = 0.5 (^k, dotted line) and a = 1 (♦, 
solid line). Again, the classical threshold is shown for com- 
parison (■). 



write the form of the input state-encoded cluster state 



|^) D k3,4,5 + &|l)l(h) D -^ 



11 |M) D )2,3,4,5, 

(24) 

where \rj) IJ (|/x) ) is the part of the noisy subcluster state 
which has the first qubit in |0) (|1)) and the qubit labels 
have been explicitly introduced. The measurement pat- 
tern in Fig. 9 (a) together with the specified assumption 
about the set of outcomes, leads directly to the final state 
present on qubit 5 

)' |Vaut) = U (l + + e iv - e ^ +u+e A + ^1 - a 2 (l + e 



a [ 1 + e l « - e ^+ 8 3) 



e«(|0)-e^|l)) 5 . 



(25) 



A. Arbitrary rotation 

The approach used to simulate arbitrary rotations 
Ur € SU{2) [1] requires the decomposition of the ro- 
tation in terms of the elements of the Euler-angle vec- 
tor fi = (C,^, C) as Ur = R X R Z R X , where the elemen- 
tary rotations around the j-axis (j — x, z) are given by 
Rj = cxp (— iOi(Tj/2). In order to simulate the action of 
Ur on the input state |^i n ), we require a five-qubit lin- 
ear cluster state and the pattern shown in Fig. 9 (a) [1]. 
The pattern consists of single-qubit measurements along 
directions in direct relation to the Euler angles. These 
relations can be seen very easily if one pictures Fig. 9 (a) 
as a concatenation of four BBB^s. In this case we have 
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1 Hir 7 - 3 
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(23) 



where a' 

ty{0) 



- ("1 



a, & 
7 with tit = 



(-lyi (3 and i 



ry(0) 



(-1) S 1 ^3 <y W itn Uj, — (7/ ~ " Ox 

The measurement bases have now become adaptive, how- 
ever this does not pose a problem as long as we mea- 
sure in the order of the qubits in the cluster. Then 
by redefining the angles for the measurement bases 
as a = (-l)«? ,co) (-£),/? = (-iy7 (a) (-v) and 7 = 

(— l) s i M< '® s 3 vW) (— W e simulate the correct rotation 
Ur. The procedure is now used to simulate the rota- 
tion of the input state \ipi n ) = a\0) +b\l), (b = — a 2 ) 
via a five-qubit noisy linear cluster. We assume the post- 
selection of the measurement results such that we may 
retain only the case corresponding to the measurement 
of the tensorial product |+)2 |+)f |+)2- It is easy to 



This should be compared to the ideally rotated state 
Ur |V>in>- The fidelity F R5 = | (^ out | U R |^ in ) | 2 , averaged 
over any possible input state and over individual Gaus- 
sian distributions associated with each unwanted 9j, is 
shown in Fig. 10 (a) for different values of the Euler 
angles (we label this average benchmark as Frs). The 
initial state dependences are washed out by assuming a 
uniform distribution for each value of a, while the Gaus- 
sians are all centered on 6j = and have equal standard 
deviation o. It is evident that the average gate fidelity 
has a considerable dependence on the particular rotation 
we want to simulate. In this case however, an average 
over all the rotation angles is meaningless, as the choice 
of the elements of the vector $7 is imposed by the spe- 
cific computation protocol desired. The general trend, 
irrespective of the Euler angles, is that the gate fidelity 
is reduced by the increase of randomness in the system 
given by the value of a. The reduction can be quite con- 
siderable even within moderate deviations from the ideal 
values 8j = 0, Vj (see the ■ or the ♦ case in Fig. 10 (a) 
for example, which suffer an average fidelity reduction of 
about 20%, for o — 0.4). 

However, this kind of analysis is unable to highlight the 
effect of redundant qubits initially present in the physical 
configuration of a cluster. These must be removed before 
the effective gate simulation is performed [1]. Such a re- 
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FIG. 9: (a): The layout for a single-qubit rotation where the 
input logical state encoded on qubit 1, is rotated by Ur and 
transfered to qubit 5 after the measurements shown, (b): A 
modified configuration with two redundant qubits (qubit 3 
and 6) which are removed via a x measurements. 
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moval of qubits, not necessary for the simulation of a par- 
ticular gate, can be achieved by measuring them in the 
o~ x or cr 2 eigenbasis. For the case of measurements in the 
cr z eigenbasis, we break any entanglement between that 
qubit and the rest of the cluster. Whereas for measure- 
ments in the a x eigenbasis, when pairs of qubits are mea- 
sured, we obtain a reduced cluster state where the qubits 
no longer affect the QIP protocol being simulated [1]. 
This removal is a very important point and the influence 
of noise on QIP protocols has only been partially clari- 
fied by our study about quantum state transfer. Indeed, 
a measurement of a qubit in a state affected by a set 
of Oj's, spreads noise throughout the cluster. More pre- 
cisely, the unwanted phase 6j which was attached to the 
j—th qubit, eliminated from the cluster by a measure- 
ment, is then inherited by the surviving qubits surround- 
ing the j-th one. This leaves us with a smaller cluster 
state which is plagued by a larger number of unwanted 
phases. This noise inheritance effect must obviously be 
kept to a minimum. 

Moreover, this effect does not depend on the particular 
angle chosen for the measurement basis and noise inher- 
itance appears after any measurement (belonging to a 
legitimate measurement pattern) on a cluster state. It 
thus becomes an intrinsic feature of the inherently noisy 
cluster state generation we are addressing [13]. A way 
to give an explicit account of the inheritance effect is by 
fictitiously modifying the configuration in Fig. 9 (a) as 
shown in panel (b), where qubits 3 and 6 must be con- 
sidered as redundant. After their elimination, via mea- 
surements in the a x -eigenbasis (which do not break the 
channels between the surviving qubits 2 & 4 and 5 & 7) , 
the physical layout is exactly the five-qubit linear clus- 
ter considered in panel (a). We have calculated the gate 
fidelity Fry after the removal of these qubits and the 
correction of the resulting cluster state via local oper- 
ations, as if a Hadamard gate has been performed be- 
tween qubits 3 & 4 and 6 & 7. Analyzing the effect of 
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FIG. 10: (a): The fidelity of one-qubit rotations decomposed 
using the Euler angles (£, v, £), when our model of noise is 
considered. The fidelity is plotted against a common standard 
deviation a on the unwanted phases for various Euler angles. 
We have considered £ = 7r/4, v = £ = (♦, solid line), £ = 
v = tt/2, i = (★, dashed line), C = 0, v = £ = vr/4 ■ 
dot-dashed line) and £ = £ = 0, v = -k (A, dotted line); 
(b): The difference Frs — Fm between the gate fidelity of 
the configurations in Figs. 9 (a) and (b). 



noise spreading through the measurements, we find that 
the average of Frj over the input states and individ- 
ual Gaussian distributions (Fry) is always smaller than 
Fr5. This is shown in Fig. 10 (b), where we have plotted 
the difference Frq — Frj (always positive) against a, for 
the rotation angles considered in Fig. 10 (a). The dif- 
ferences between the two cluster configurations can eas- 
ily exceed 10% and are larger for increased randomness 
in the cluster. Their behaviors as the Euler angles are 
changed, are almost uniform until a ~ 0.6. Then for 
increased randomness, the specific way in which the in- 
herited phases are distributed within the structure of the 
cluster state becomes relevant and some discrepancies oc- 
cur. The message here is that there is a counterintuitive 
dependence of the gate fidelity on the specific angles of 
rotation: some rotations are more exposed to our model 
of noise than others. 



B. CNOT gate 

In order to complete our analysis of QIP protocols in 
the QC C scheme with intrinsic noise, we address the sim- 
ulation of a CNOT as the paradigmatic example of an 
entangling two-qubit gate [3]. We first examine the be- 
havior of the original proposal [1] under noisy conditions 
and in the presence of a redundant qubit. Then, we will 
compare it with other configurations for CNOT simula- 
tion, where the resource requirements can be dramati- 
cally reduced down to no more than four qubits [13]. 

The input states will be denoted \Qi) = a|0) + 6|1) 
and IQ2) = c|0) + d|l) which encode the control and 
target state respectively. We consider the fifteen-qubit 
bidimensional cluster state (the squashed-l configuration) 
whose layout and measurement pattern are shown in 
Fig. (11) (a). Its workings in terms of stabilizer for- 
malism can be found in [1]. Here, we are interested in 
the concatenation technique, which gives an immediate 
picture of the equivalent quantum gates simulated by this 
cluster configuration. Using the BBB's and EBB's in- 
troduced in Section III, it is straightforward to derive 
the equivalent quantum circuit in Fig. 11 (b). The role 
played by the EBB introduced in Fig. 6, bridging the 
two otherwise independent subclusters, is crucial here. 
The byproduct operator can be evaluated starting from 
this equivalent decomposition and explicitly found to be 
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where sf 

We now consider the simulation of a CNOT with a noisy 
squashed-I cluster. The technique already exploited in 
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FIG. 11: (a): Squashed-I configuration for a CNOT simula- 
tion. The input (output) control and target logical qubit are 
qubits 1 (7) and 9 (15) respectively, with qubit 8 as a bridging 
qubit. (b): The equivalent quantum circuit as a concatena- 
tion of BBB's and EBB's. 
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Explicit form in the computational basis 
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TABLE I: The notation used in Eq. (28) for the noisy 
squashed-I cluster state used to simulate a CNOT. 



Section VI A in order to write an implicit (but manage- 
able) expression for the global state of a cluster in a given 
configuration is used again and the starting point for our 
calculations is the state 



I+)8®IV)^ + I2®b)f 3 ^15- ( 27 ) 



\r/) D + \fi) D , where \ip) D is the 



l<P)f-4® l^)^->7 ' 

Here \ip) D = |^) D + |x) J 
part of I ip) D which has the last qubit in |0), while \x) D is 
the part with the last qubit in |1). The explicit form of 

W) i \x) D ! \v) D an( i \^) D is given in Tab. I. In Eq. (27) 
the encoded logical input state of the control and target 
qubits, as well as the state of the bridging qubit 8, have 
been properly singled out. The 1 — > 4 (9 — > 12) symbol 
means that qubits 1, 2, 3, 4 (9, 10, 11, 12) are involved. By 



applying the entangling operations S^ 8 , S^ lz , S*^ and 
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Sp' 1 , the state of the cluster becomes 



5^7 

D 



WL4\^7 + \X)^(\V) D - 



|0) 8 ® 



13->15 



13->15 

(28) 



Here, Of (6%) denotes the unwanted phase introduced 
when the entanglement between qubits 4 and 5 (qubits 
4 and 8) is considered. Using this notation, the entan- 
gling operation which glues together two subclusters is 



considerably simplified. For instance, gluing |y)^_ >4 and 
1^)5^7 leads to 



\xf)l^4® 



|x)l-4 \lA 



5^7 ' 



5-+7 

(29) 



which is a seven-qubit linear cluster state. It should be 
noted that this expression has been obtained without ac- 
tually writing the subcluster states in the computational 
basis. This method is therefore space-saving and compu- 
tationally useful. On the other hand, by dealing each 
time with small subcluster states, it is handy to find 
out the explicit form of the state of the output logical 
qubits after the application of the appropriate measure- 
ment pattern. It is easy to build up a table for the trans- 
formations occurring after the completed measurements. 

Following these lines, the gate fidelity i^cNOT can be 
explicitly evaluated by hand. A plot is given in Fig. 12 
(a) against both 6j — 9 (Vj involved in Eq. (28)) and 
a = c (for convenience) with the normalization condi- 
?, d 



tions b = VT 



vT 



As in the case of single- 



qubit rotations, we have considered the post-selection of 
the event corresponding to s\ — = among the set 
of outcomes resulting from the measurements. In this 
case, the decoding operator is U^({0}) = a^P ® l^ 15 -*. 
Behaviors qualitatively similar to Fig. 12 (a) can be ob- 
served for any other choice of the relation between a and 



Gate Fidelity 




()?2.0:4-().6 0.8 1 



FIG. 12: (a): Fidelity for a squashed-I CNOT plotted against 
the unwanted phase 9 and the input control-state coefficient a. 
In this plot a = c. (b) : Gate fidelities for different simulations 
of a CNOT. From top to bottom curve, we show the fidelity 
of the four-qubit CNOT of Fig. 5, the helix configuration, the 
squashed-I and the squashed-I with an additional bridging 
qubit. 
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c with this plot having the merit of showing an almost 
uniform behavior of -Fcnot against a, for a fixed 9 and a 
fast decay of the gate fidelity is found for non-zero values 
of the unwanted phases. Near 9 — 0.6 and a = c = 0.5, 
-Fcnot — 0.2 is found and is never exceeded for fixed 
a, whatever the choice for the relation between a and 
c. By showing the fidelity behavior in Fig. 12 (a), we 
produce a significant example of the performances of the 
squashed-I CNOT simulation in the presence of our model 
of noise, showing that the one-way model has to face a 
considerable decay in the two-qubit gate fidelity. These 
conclusions are strengthened by the calculation of the 
average fidelity over Gaussian distributions for the un- 
wanted phases shown in the dot-dashed line (♦ symbol) 
in Fig. 12 (b). Although a more complete analysis re- 
quires an average over all the outcome configurations [23] , 
this example provides sufficient physical insight. 

We can study the way in which the noise inheritance 
attacks the CNOT fidelity by analyzing the simple but 
interesting example of a modified squashed-I where the 
bridging zone between the control and target subclus- 
ters includes an additional redundant qubit. We remove 
this qubit from the cluster (via a measurement in the 
tTx-eigenbasis) retrieving the fifteen-qubit cluster state 
through appropriate local operations on qubit 12. Fol- 
lowing the lines depicted above i.e. attaching Gaussian 
distributions of standard deviation a to the unwanted 
phases appearing in the noisy cluster state, we can ver- 
ify that the average gate fidelity is spoiled, as shown in 
Fig. 12 (b) (A, dashed line) from the expected spread of 
additional phases relative to the removed qubit [13]. 



C. Alternative routes to CNOT 

Our discussion about noise inheritance reinforces the 
view we introduced previously, concerning the impor- 
tance of keeping the number of qubits in a cluster as 
low as possible. The analysis of the squashed-I layout re- 
vealed it to be rather prone to the effect of the intrinsic 
noise model. We ascribe this to the expensive nature of 
the configuration in terms of the number of qubits in the 
cluster. Thus there is a necessity for looking at differ- 
ent ways in which an entangling two-qubit gate can be 
simulated through cluster states. Here, we would like to 
provide an example of alternative cluster configurations 
which are able to simulate a CNOT gate involving less 
qubits than the squashed-I. 

We consider the qubit layout and measurement pat- 
tern sketched in Fig. 13 (a), which simulates a CNOT 
using ten qubits. We will refer to this configuration as 
the helix configuration. The scheme is based on the sim- 
ulation of a CrPHASE gate [3] (within the dashed box) 
and realizes the transformation H(C7rPHASE)H = CNOT. 
A more detailed analysis of the gate simulation, using the 
concatenation technique is given in Fig. 13 (b), where a 
relabelling of the logical output qubits, effectively equiv- 
alent to a SWAP gate is required. The crucial feature is 
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FIG. 13: (a): Helix layout and measurement pattern for a 
CNOT simulated through a CRPHASE gate (within the dashed 
square) and two Hadamard gates involving qubits 6 & 9 and 
8 & 10. (b): Equivalent quantum circuit drawn by exploiting 
the concatenation of eight BBB\ 's and one box-shaped EBB 
involving qubits 1,2,3,4 (see Fig. 4). 



the use of a box cluster (see Fig. 4) . 

It is possible to carry out an analytic calculation of the 
dynamics of a noisy helix cluster state using the same 
technique highlighted previously. The resulting average 
gate fidelity is shown in Fig. 12 (b) (■, dotted line). 
As before, we assume the post-selection of the event cor- 
responding to the set of outcomes sf = (i.e. all the 
measured qubits are found in [+)). This results in no 
decoding operator at the end of the procedure, a clear 
advantage with respect to the squashed-I CNOT. The 
usual individual Gaussians have been considered, with 
tj as their standard deviation and a noticeable improve- 
ment in the gate fidelity is observed, compared to the 
squashed-I configuration. 

The situation can be further improved by looking at 
the squashed-I cluster configuration and examining more 
closely the simulation performed there. It is straightfor- 
ward to recognize that the non-local nature of the CNOT 
gate in Fig. 11 (b) is all in the C^ PHASE gate sandwiched 
by the Hadamards on the IQ2) — IQ2) (target) line. The 
remainder of the circuit realizes local operations on the 
control qubit, which are unnecessary for the CNOT simu- 
lation. Thus, there is a consistent redundancy in this clus- 
ter configuration. Stripping the squashed-I cluster bare 
to eliminate the unnecessary local parts, leads directly to 
the EBB already introduced in Fig. 6 at its very core. 
But even this is unnecessary as we already know that it 
is possible to obtain a better configuration which pro- 
vides an even more economical configuration for a CNOT 
simulation. This is the four-qubit EBB discussed in Sec- 
tion III Fig. 5, which naturally simulates a CNOT (with 
output states in the cr x -eigenbasis). Our aim is therefore 
confronting the performance of the noisy version of this 
simple configuration with the other simulations treated 
so far. The results are shown in Fig. 12 (b) (-k and solid 
line). Evidently the gate fidelity corresponding to this 
linear layout is vastly superior to any other case treated 
so far [13]. This four-qubit EBB CNOT is important for 
two essential reasons. First, it can be seen as the ulti- 
mate confirmation that in a noisy scenario, those features 
which the ideal one-way model takes for granted, i.e. the 
management of an arbitrarily large cluster and the in- 
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nocuousness of the measurements performed in order to 
process the encoded information, sensibly affect a com- 
putational task when imperfections are embedded in the 
cluster state. Second, it is evident that more economical 
configurations for gate simulation are required if reliable 
computation is to be performed. The four-qubit linear 
EBB actually provides an economical CNOT simulation 
in which the number of parameters involved corresponds 
to the number of qubits in the cluster resource. These 
points, already extensively commented in [13], have been 
enriched here by a more detailed analysis. 

VII. REMARKS 

We have addressed several very important points in the 
cluster state model for QC. The construction and analy- 
sis of cluster configurations suitable for the simulation of 
desired quantum circuits, has been simplified by the in- 
troduction of a class of elementary BBB's. We have seen 
how they can be concatenated together in order to model 
even complicated circuit configurations. Our strategy al- 
lows one to immediately figure out the form of the de- 
coding operator in order to eliminate the effects of the 
measurement-induced randomness in a simulation. We 



believe this represents a valuable tool for the purposes 
of configuration-design, which is helpful both theoreti- 
cally and experimentally. Indeed, by properly designing 
the cluster state resource, it is possible to minimize the 
number of redundant qubits in a circuit simulation. This 
point is related to the second task accomplished by this 
paper, namely the detailed analysis of QIP using intrin- 
sic imperfections in the generation of cluster configura- 
tions. We have shown that uncontrollable randomness in 
the qubit-qubit interactions which create a cluster state, 
affect both communication and computation protocols 
based on the one-way model. A direct consequence of 
our analysis is that in the processing of information en- 
coded in a cluster state, both the number of qubits in- 
volved and the measurements to be performed must be 
carefully managed. Our study paves the way toward re- 
search of gate simulations performed by using only small 
clusters of just a few qubits. The existence of such eco- 
nomic configurations has been further commented here. 
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